A comprehensive, multi-process box-model approach to glacial-interglacial carbon cycling

The canonical question of which physical, chemical or biological mechanisms were responsible for oceanic uptake of atmospheric CO2 during the last glacial is yet unanswered. Insight from paleo proxies has led to a multitude of hypotheses but none so far have been convincingly supported in three dimensional numerical modelling experiments. The processes that influence the CO2 uptake and export production are inter-related and too complex to solve conceptually while complex numerical models are time consuming and expensive to run which severely limits the combinations of mechanisms that can be explored. Instead, an intermediate inverse box model approach is used here in which the whole parameter space is explored. The glacial circulation and biological production states are derived from these using proxies of glacial export production and the need to draw down CO2 into the ocean. We find that circulation patterns which explain glacial observations include reduced Antarctic Bottom Water formation and high latitude mixing and to a lesser extent reduced equatorial upwelling. The proposed mechanism of CO2 uptake by an increase of eddies in the Southern Ocean, leading to a reduced residual circulation, is not supported. Regarding biological mechanisms, an increase in the nutrient utilization in either the equatorial regions or the northern polar latitudes can reduce atmospheric CO2 and satisfy proxies of glacial export production. Consistent with previous studies, CO2 is drawn down more easily through increased productivity in the Antarctic region than the sub-Antarctic, but that violates observations of lower export production there

Interhemispheric coupling, the West Antarctic Ice Sheet and warm Antarctic interglacials

Ice core evidence indicates that even though atmospheric CO2 concentrations did not exceed 300 ppm at any point during the last 800 000 years, East Antarctica was at least 3–4 C warmer than preindustrial (CO2 280 ppm) in each of the last four interglacials. During the previous three interglacials, this anomalous warming was short lived (3000 years) and apparently occurred before the completion of Northern Hemisphere deglaciation. Hereafter, we refer to these periods as “Warmer than Present Transients” (WPTs). We present a series of experiments to investigate the impact of deglacial meltwater on the Atlantic Meridional Overturning Circulation (AMOC) and Antarctic temperature. It is well known that a slowed AMOC would increase southern sea surface temperature (SST) through the bipolar seesaw and observational data suggests that the AMOC remained weak throughout the terminations precedingWPTs, strengthening rapidly at a time which coincides closely with peak Antarctic temperature. We present two 800 kyr transient simulations using the Intermediate Complexity model GENIE-1 which demonstrate that meltwater forcing generates transient southern warming that is consistent with the timing of WPTs, but is not sufficient (in this single parameterisation) to reproduce the magnitude of observed warmth. In order to investigate model and boundary condition uncertainty, we present three ensembles of transient GENIE-1 simulations across Termination II (135 000 to 124 000 BP) and three snapshot HadCM3 simulations at 130 000 BP. Only with consideration of the possible feedback of West Antarctic Ice Sheet (WAIS) retreat does it become possible to simulate the magnitude of observed warming

A model-based constraint on CO2 fertilisation

We derive a constraint on the strength of CO2 fertilisation of the terrestrial biosphere through a "top-down" approach, calibrating Earth system model parameters constrained by the post-industrial increase of atmospheric CO2 concentration. We derive a probabilistic prediction for the globally averaged strength of CO2 fertilisation in nature, for the period 1850 to 2000 AD, implicitly net of other limiting factors such as nutrient availability. The approach yields an estimate that is independent of CO2 enrichment experiments. To achieve this, an essential requirement was the incorporation of a land use change (LUC) scheme into the GENIE Earth system model. Using output from a 671-member ensemble of transient GENIE simulations, we build an emulator of the change in atmospheric CO2 concentration change since the preindustrial period. We use this emulator to sample the 28-dimensional input parameter space. A Bayesian calibration of the emulator output suggests that the increase in gross primary productivity (GPP) in response to a doubling of CO2 from preindustrial values is very likely (90% confidence) to exceed 20%, with a most likely value of 40–60%. It is important to note that we do not represent all of the possible contributing mechanisms to the terrestrial sink. The missing processes are subsumed into our calibration of CO2 fertilisation, which therefore represents the combined effect of CO2 fertilisation and additional missing processes. If the missing processes are a net sink then our estimate represents an upper bound. We derive calibrated estimates of carbon fluxes that are consistent with existing estimates. The present-day land–atmosphere flux (1990–2000) is estimated at −0.7 GTC yr−1 (likely, 66% confidence, in the range 0.4 to −1.7 GTC yr−1). The present-day ocean–atmosphere flux (1990–2000) is estimated to be −2.3 GTC yr−1 (likely in the range −1.8 to −2.7 GTC yr−1). We estimate cumulative net land emissions over the post-industrial period (land use change emissions net of the CO2 fertilisation and climate sinks) to be 66 GTC, likely to lie in the range 0 to 128 GTC

Impact of negative and positive CO2 emissions on global warming metrics using an ensemble of Earth system model simulations

The benefits of implementing negative emission technologies in the global warming response to cumulative carbon emissions until the year 2420 are assessed following the shared socioeconomic pathway (SSP) 1-2.6, the sustainable development scenario, with a comprehensive set of intermediate-complexity Earth system model integrations. Model integrations include 86 different model realisations covering a wide range of plausible climate states. The global warming response is assessed in terms of two key climate metrics: the effective transient climate response to cumulative CO2 emissions (eTCRE), measuring the surface warming response to cumulative carbon emissions and associated non-CO2 forcing, and the effective zero emissions commitment (eZEC), measuring the extent of any continued warming after net-zero CO2 emissions are reached. The transient climate response to cumulative CO2 emissions (TCRE) is estimated as 2.2 K EgC−1 (median value) with a 10 %–90 % range of 1.75 to 3.13 K EgC−1 in 2100, approximated from the eTCRE by removing the contribution of non-CO2 forcing. During the positive emission phase, the eTCRE decreases from 2.71 (2.0 to 3.65) to 2.61 (1.91 to 3.62) K EgC−1 due to a weakening in the dependence of radiative forcing on atmospheric carbon, which is partly opposed by an increasing fraction of the radiative forcing warming the surface as the ocean stratifies. During the net negative and zero emission phases, a progressive reduction in the eTCRE to 2.0 (1.39 to 2.96) K EgC−1 is driven by the reducing airborne fraction as atmospheric CO2 is drawn down mainly by the ocean. The model uncertainty in the slopes of warming versus cumulative CO2 emissions varies from being controlled by the radiative feedback parameter during positive emissions to being affected by carbon-cycle parameters during net negative emissions, consistent with the drivers of uncertainty diagnosed from the coefficient of variation of the contributions in the eTCRE framework. The continued warming after CO2 emissions cease and remain at zero gives a model mean eZEC of −0.03 K after 25 years, which decreases in time to −0.21 K at 90 years after emissions cease. However, there is a spread in the ensemble with a temperature overshoot occurring in 20 % of the ensemble members at 25 years after cessation of emissions. If net negative emissions are included, there is a reduction in atmospheric CO2 and there is a decrease in temperature overshoot so that the eZEC is positive in only 5 % of the ensemble members. Hence, incorporating negative emissions enhances the ability to meet climate targets and avoid risk of continued warming after net zero is reached

Stability of the trapped nonconservative Gross-Pitaevskii equation with attractive two-body interaction

The dynamics of a nonconservative Gross-Pitaevskii equation for trapped atomic systems with attractive two-body interaction is numerically investigated, considering wide variations of the nonconservative parameters, related to atomic feeding and dissipation. We study the possible limitations of the mean field description for an atomic condensate with attractive two-body interaction, by defining the parameter regions where stable or unstable formation can be found. The present study is useful and timely considering the possibility of large variations of attractive two-body scattering lengths, which may be feasible in recent experiments.Comment: 6 pages, 5 figures, submitted to Physical Review

Coupling climate and economic models in a cost-benefit framework: a convex optimization approach

In this paper we present a general method, based on a convex optimisation technique, that facilitates the coupling of climate and economic models in a cost-benefit framework. As a demonstration of the method, we couple an economic growth model à la Ramsey adapted from DICE-99 with an efficient intermediate complexity climate model, C-GOLDSTEIN, which has highly simplified physics, but fully 3-D ocean dynamics. As in DICE-99 we assume that an economic cost is associated with global temperature change: this change is obtained from the climate model which is driven by the GHG concentrations computed from the economic growth path. The work extends a previous paper in which these models were coupled in cost-effectiveness mode. Here we consider the more intricate cost-benefit coupling in which the climate impact is not fixed a priori. We implement the coupled model using an oracle-based optimisation technique. Each model is contained in an oracle which supplies model output and information on its sensitivity to a master program. The algorithm Proximal-ACCPM guarantees the convergence of the procedure under sufficient convexity assumptions. Our results demonstrate the possibility of a consistent, cost-benefit, climate-damage optimisation analysis with a 3-D climate model

Height and risk of death among men and women: aetiological implications of associations with cardiorespiratory disease and cancer mortality

OBJECTIVES: Height is inversely associated with cardiovascular disease mortality risk and has shown variable associations with cancer incidence and mortality. The interpretation of findings from previous studies has been constrained by data limitations. Associations between height and specific causes of death were investigated in a large general population cohort of men and women from the West of Scotland. DESIGN: Prospective observational study. SETTING: Renfrew and Paisley, in the West of Scotland. SUBJECTS: 7052 men and 8354 women aged 45-64 were recruited into a study in Renfrew and Paisley, in the West of Scotland, between 1972 and 1976. Detailed assessments of cardiovascular disease risk factors, morbidity and socioeconomic circumstances were made at baseline. MAIN OUTCOME MEASURES: Deaths during 20 years of follow up classified into specific causes. RESULTS: Over the follow up period 3347 men and 2638 women died. Height is inversely associated with all cause, coronary heart disease, stroke, and respiratory disease mortality among men and women. Adjustment for socioeconomic position and cardiovascular risk factors had little influence on these associations. Height is strongly associated with forced expiratory volume in one second (FEV1) and adjustment for FEV1 considerably attenuated the association between height and cardiorespiratory mortality. Smoking related cancer mortality is not associated with height. The risk of deaths from cancer unrelated to smoking tended to increase with height, particularly for haematopoietic, colorectal and prostate cancers. Stomach cancer mortality was inversely associated with height. Adjustment for socioeconomic position had little influence on these associations. CONCLUSION: Height serves partly as an indicator of socioeconomic circumstances and nutritional status in childhood and this may underlie the inverse associations between height and adulthood cardiorespiratory mortality. Much of the association between height and cardiorespiratory mortality was accounted for by lung function, which is also partly determined by exposures acting in childhood. The inverse association between height and stomach cancer mortality probably reflects Helicobacter pylori infection in childhood resulting inor being associated withshorter height. The positive associations between height and several cancers unrelated to smoking could reflect the influence of calorie intake during childhood on the risk of these cancers

System complexity and policy integration challenges: The Brazilian Energy- Water-Food Nexus

The Energy-Water-Food Nexus is one of the most complex sustainability challenges faced by the world. This is particularly true in Brazil, where insufficiently understood interactions within the Nexus are contributing to large-scale deforestation and land-use change, water and energy scarcity, and increased vulnerability to climate change. The reason is a combination of global environmental change and global economic change, putting unprecedented pressures on the Brazilian environment and ecosystems. In this paper, we identify and discuss the main Nexus challenges faced by Brazil across sectors (e.g. energy, agriculture, water) and scales (e.g. federal, state, municipal). We use four case studies to explore all nodes of the Nexus. For each, we analyse data from economic and biophysical modelling sources in combination with an overview of the legislative and policy landscape, in order to identify governance shortcomings in the context of growing challenges. We analyse the complex interdependence of developments at the global and local (Brazilian) levels, highlighting the impact of global environmental and economic change on Brazil and, conversely, that of developments in Brazil for other countries and the world. We conclude that there is a need to adjust the scientific approach to these challenges as an enabling condition for stronger science-policy bridges for sustainability policy-making. © 2019 The Author(s

Oracle-based optimization applied to climate model calibration

In this paper, we show how oracle-based optimization can be effectively used for the calibration of an intermediate complexity climate model. In a fully developed example, we estimate the 12 principal parameters of the C-GOLDSTEIN climate model by using an oracle- based optimization tool, Proximal-ACCPM. The oracle is a procedure that finds, for each query point, a value for the goodness-of-fit function and an evaluation of its gradient. The difficulty in the model calibration problem stems from the need to undertake costly calculations for each simulation and also from the fact that the error function used to assess the goodness-of-fit is not convex. The method converges to a Fbest fit_ estimate over 10 times faster than a comparable test using the ensemble Kalman filter. The approach is simple to implement and potentially useful in calibrating computationally demanding models based on temporal integration (simulation), for which functional derivative information is not readily available

Producing Policy-relevant Science by Enhancing Robustness and Model Integration for the Assessment of Global Environmental Change

We use the flexible model coupling technology known as the bespoke framework generator to link established existing modules representing dynamics in the global economy (GEMINI_E3), the energy system (TIAM-WORLD), the global and regional climate system (MAGICC6, PLASIM-ENTS and ClimGEN), the agricultural system, the hydrological system and ecosystems (LPJmL), together in a single integrated assessment modelling (IAM) framework, building on the pre-existing framework of the Community Integrated Assessment System. Next, we demonstrate the application of the framework to produce policy-relevant scientific information. We use it to show that when using carbon price mechanisms to induce a transition from a high-carbon to a low-carbon economy, prices can be minimised if policy action is taken early, if burden sharing regimes are used, and if agriculture is intensified. Some of the coupled models have been made available for use at a secure and user-friendly web portal. © 2018 The Author